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■ Abstract 

(N 

In this paper, we obtain the ladder operators and associated compatibility con- 

■ ditions for the type I and the type II multiple orthogonal polynomials. These 

■ ladder equations extend known results for orthogonal polynomials and can be used 
to derive the differential equations satisfied by multiple orthogonal polynomials. 
Our approach is based on Riemann-Hilbert problems and the Christoffel-Darboux 
formula for multiple orthogonal polynomials, and the nearest-neighbor recurrence 
relations. As an illustration, we give several explicit examples involving multiple 
Hermite and Laguerre polynomials, and multiple orthogonal polynomials with ex- 
ponential weights and cubic potentials. 
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1 Introduction 



1.1 Multiple orthogonal polynomials 

^ \ Multiple orthogonal polynomials are polynomials of one variable which are defined by 

orthogonality relations with respect to r different measures jix, \x-i-, ■ ■ ■ , where r > 1 
is a positive integer. As a generalization of orthogonal polynomials [T5l [25] [Mj . multiple 
orthogonal polynomials originated from Hermite-Pade approximation in the context of 
irrationality and transcendence proofs in number theory. They were further developed in 
approximation theory, cf . j2j HJ HH1 EEJ [33] and surveys [HI [35J [36] . During the past few 
years, multiple orthogonal polynomials have also arisen in a natural way in certain models 
from mathematical physics, including random matrix theory, non- intersecting paths, etc.; 
we refer to [301 EI] and references therein for the progress of this subject. 

Let n = (ni, n 2 , . . . , n r ) G N r be a multi-index of size |n| = n 1 +n 2 +. . .+n r and suppose 
fix, A*2, ■ • • , A*r are r measures with supports on certain simple curves in the complex plane. 
There are two types of multiple orthogonal polynomials. The type I multiple orthogonal 
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polynomials are given by the vector (A^i, . . . , where Aftj is a polynomial of degree 

< rij — 1, for which 



We use the normalization 



y x fe y4^(x) djt/j(a;) = 0, fc = 0, 1, 



Inl - 2. 



x |n| Afjj(x) dfij(x) — 1. 



(1-1) 



;i.2) 



The equations (ll.ip — (ll.2p are a linear system of |n| equations for the unknown coefficients 
of Afj i, . . . , Ati^ r . This system has a unique solution if the matrix of mixed moments 
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where 



x n dfj,j(x) 



J 



is not singular, in which case we call the multi-index n a normal index. The type II 
multiple orthogonal polynomial is the monic polynomial Ph(x) = x' n ' + • • • of degree 
for which 



n 



J P H (x)x k d^i(x) = 0, k = 0, 1, ... 



;i.3) 



Pft(x)x k dji r {x) = 0, 



= 0, 1, . . . , n r — 1. 



This gives a linear system of |n| equations for the |n| unknown coefficients of Pa and the 
matrix of this linear system is the transpose Mj of (11. ip . hence P^ exists and is unique 
whenever n is a normal index. 

Suppose that the r measures //i, . . . , /i r are all absolutely continuous with respect to 
a measure /x and that dfij(x) = Wj(x) dfi(x). The type I and II multiple orthogonal 
polynomials satisfy the following biorthogonality: 



Pn{x)Qm( X ) d K X ) 



if in < n, 

if \n\ < \m\ — 2, 



:i-4) 



if \m\ 



\n\ + 1, 



where 

r 

Qn{x) = } j Aft,j(x)w j (x), (1.5) 

3=1 

see [281 Theorem 23.1.6]. 

For more information about multiple orthogonal polynomials, we refer to Aptekarev et 
al. [31 Coussement and Van Assche [38J, Nikishin and Sorokin (331 Chapter 4, §3], and 
Ismail [2S1 Chapter 23]. Throughout this paper, we shall assume that all multi- indices 
are normal. 
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1.2 Recurrence relations 

There are several recurrence relations for the type I and the type II multiple orthogonal 
polynomials. For type II, we have the following nearest-neighbor recurrence relations 

r 

xPn(x) = Pfi-^g 1 {x) + bfi t iPn(x) + ^ dftjPft-g. {x) , 

3=1 

: (1.6) 

r 

xPfi{x^ P-tl'Sr (•^') "I - bfi^rPn^X^) -\- ^ ^ ^n,j Pft—ej (»e) ; 

where = (0, . . . , 0, 1, 0, . . . , 0) is the j-th standard unit vector with 1 on the j-th en- 
try, (a^i, . . . , a^ r ) and (fe^i, . . . , b^ r ) are the recurrence coefficients. These recurrence 
relations were derived in [28| Theorem 23.1.11] where it was shown that 



x nj Pft(x) djjj(x) 
/ x ni ~ 1 P a -s j {x)dii j (x) 



and 

bftj= I xPft(x)Q n+ g (x) dfi(x), (li 



for j = 1, . . . , r, where Q,? is given in (11.51) . 

We have similar recurrence relations for the type I multiple orthogonal polynomials: 

r 

Qri— e*i "I - bfl—g 1 \Qfl{x) -\- ^ ^ Q"n,jQn+ej (^Q ; 

i=i 

: (i.9) 

r 
3=1 

Observe that the same recurrence coefficients (a^i, . . . , a^ jr ) are used but that there is a 
shift in the other recurrence coefficients (bfts lt i, • • • , &n-e r ,r)- 

If we multiply (11.61) by Qn(x) and integrate, then the biorthogonality (II. 4p gives 

r 

xPji{x^C^ f[(yX^ d^i{x^ ^ ^ ^ft,j ' 

3=1 

The orthogonality properties of P^ in (11. 3B imply that 

Pfi(a:)Qn+e,-(^) dfi(x) = / P s (a;)^i+aj(^) dfij(x), 



hence, by (11.41) . 

l = K H+ g jtj I x n iP H {x)dnj{x), (1.10) 



where Aa t j(x) = K^x n ' 1 + • • • . This, together with (jl.7p . provides an alternative repre- 
sentation of a f , 



J x nj Pfi(x) dfi 



! 1 ' ' ~ r ' J \ *^ / ^'Tl J / -j -j -j \ 

J X n i ^Pfi—g. [x) dflj(x) Hfi+ejj 

1.3 Ladder equations for orthogonal polynomials 

Given a single positive measure w(x)dx, it is well-known that there exist monic orthogonal 
polynomials P n {x) of degree n in x such that 



P m {x)P n (x)w(x) dx = h n 5 m)n , h n >0, m,n = 0,1,2, — (1-12) 

They are characterized by the three-term recurrence relation: 

xP n (x) = P n+1 {x) + a n P n {x) + /3„P„_i(x), (1.13) 

where 

a n = — I xP^{x)w{x) dx, (3 n = ( xP n {x)P n ^\{x)w{x) dx, (1-14) 

and the initial condition is taken to be (3qP-i{x) := 0. Suppose that w vanishes at the 
end points of the orthogonality interval, then it was shown in [TT] that P n satisfy the 
following ladder equations: 

^ + B n (x^J P n {x) = /3 n A n (x)P n _!(a;), (1.15) 

(J^ - B n (x) - y'(x^j P n ^(x) = -A n _ l {x)P n {x) (1.16) 

with v(a;) := — lnw(x) and 

An {x) : = / V/(X) ~ V ' iy) [P n (y)] 2 w(y)dy, (1.17) 
"n J % y 

B n (x) := J- / ^^"^ P n - 1 (y)P n (y)w(y)dy. (1.18) 

Furthermore, the functions v4 n and B n defined by (I1.17P and (II. 18ft satisfy 

B n+1 (z) + B n (z) = (z- a n )A n (z) -v'(z), (Si) 

1 + (z - a n )[B n+l (z) - B n (z)} = n+1 A n+ i(z) - /3„A„„i(z). (S 2 ) 

The conditions Si and S2 are usually called the compatibility conditions for the ladder 
equations. Similar relations were also achieved for discrete orthogonal polynomials [29] . 
g-orthogonal polynomials [12] and matrix orthogonal polynomials [201 [27]. 

The motivation of deriving such ladder equations is two-fold. On the one hand, the 
differential recurrence relations f II . 1 5 j) and (I1.16p . together with recurrence relation (I1.13p . 
will lead to a linear differential equation satisfied by P n . On the other hand, ladder 
operators has been successfully applied to establish the connections between Painleve 



4 



equations and recurrence coefficients of certain orthogonal polynomials; cf. j5J EJ [T3J EJ 
US I2D 1221 123] for recent applications. 

It is the aim of this paper to find and explore ladder relations for multiple orthogonal 
polynomials. We hope our results will be helpful in the further study of recurrence coef- 
ficients of multiple orthogonal polynomials and its relation to certain integrable systems. 
In the literature, the first attempt in this direction was done by Coussement and Van 
Assche [161 [38], where they derived raising and lowering operators for type II multiple 
orthogonal polynomials with respect to classical weights. The results there, however, are 
different from ours, since they consider relations between multiple orthogonal polynomials 
associated with different weights in the same class. 

The results of this paper will be presented in the next section. 

2 Statement of the results 

2.1 Ladder equations for multiple orthogonal polynomials 

Our first result concerns the ladder equations for the type II multiple orthogonal polyno- 
mials. 

Theorem 2.1. Let /ii, . . . ,/i r be r measures that are absolutely continuous with weights 
Wi, . . . , w r and each Wi vanishes at the endpoints of the support of fii. Suppose that all 
the indices n = (jii, . . . , n r ) G N r are normal and the functions 

{Wi, XWi, . . . , X ni_1 W7i, U>2, XW2, • • • , £ n2_1 U>2, • • • , W r , XW r , . . . , X nr ~ l W r } 

are linearly independent, then we have the following lowering equation for the type II 
multiple orthogonal polynomials: 

P' H {x) = P 9 (x) [ Pfi{t)Y^A n>k {t) V ' k{t) ~ V ' k{x) w k (t) dt 
J k=i x 

-J2a H ,P H ^(x) j P n (t) A n+gj 4t) V ' k{t) ~_ V 'f x) w k {t) dt, (2.1) 

j=l J k=l °° 

where v k (x) := — \ia.Wk(x) and a^j are the recurrence coefficients given in (II. 7p . 

The following raising equations for the type II multiple orthogonal polynomials hold: 
for i = 1 . . . , r one has 

PL*(x) = Pn{x) f Pn- Si {t)j^A^ k {t) V ' k{t) ~_ V ' k(x) w k (t) dt 
J k=i x 

/ p *sM ^^ +g ., fc (t) ^~ V 'f x) w k {t) dt 

j=l J k=l X 

-^(x)^-)p^(x), (2.2) 

where 8ij is the Kronecker delta. 
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Although our assumption on the weight functions may seem a little bit involved, they 
include the so-called Angelesco systems and AT systems that are normal for all indices; 
cf. [28| Sections 23.1.1-23.1.2] for an introduction of these two weight systems. Moreover, 
in the case of r = 1, suppose d/j,i(x) = w(x)dx, we have, in the notations of multiple 
orthogonal polynomials introduced in Sections 11.11 and 11.21 

P n _ x [x) 

P H {x) = P n (x), A HA (x) = — , bft,i = a n , an,i = fin, 

where P n , h n , a n and (3 n are given in (11.12j) and (11. 14ft . Inserting the above representations 
into (12. ip and (I2.2p . we recover equations (I1.15P and f II . 1 6 j) by noting the fact that /3 n = 
h n /h n -x. 

The results in Theorem 12.11 can be summarized in vector form. To this end, we 
introduce an (r + 1) x (r + 1) matrix 

N(x) := N{ft; x) = (JV«(x)) < y < r , (2.3) 

where 



Nij{x) 



J j i 



v' k (t)-v' k (x) 



x — t 

k=l 

r 



w k (t)dt, j = 0, 



J k=i 1 x 

(2.4) 

Here, it is understood that e*o = 0. We then have from (12.11) and (12.21) that 

P' H (x) = N(x)P R (x), (2.5) 

where 

P n (x) := (P n (x),P^ gl (x),...,Pn-g r (x)) T . (2.6) 
We have similar results for the type I multiple orthogonal polynomials. 

Theorem 2.2. Let . . . , fi r ber measures as given in Theorem \2.1\ For I — 1, ... ,r, the 

associated type I multiple orthogonal polynomials satisfy the following raising equations: 

A' ni {x) = -A Hil {x) [ P H (t)TA Htk (t) V ' kit) ~ V } ix) w k (t)dt 
J fi x ~ t 

+ J2a n>j A n+s . tl (x) [ P H _ s .(t)J2 A nAt) ^} ~_ V ' k{x) Wk (t) dt, (2.7) 

j=l J k=l 

where v k (x) := — hiw k (x) and a^j are the recurrence coefficients given in (JTT7J). 

The lowering equations for the type I multiple orthogonal polynomials are given by 

~ A nM [ Pn(t)j2 A ^(t) — ~_ V ' k(x) Mt) dt 

J k=l 
3=1 J k=l X 

-vKx^An+g^x), (2.8) 

for i — 1, . . . , r. 
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In vector form, Theorem 12.21 reads 



d_ 

dx 



( -A H j(x) \ 

a n,l^-n+e 1 ,l{x) 
\ a n,rA-n+e r ,l\ x ) J 



N (x) 



( Mi( x ) \ 

~ a n,l^-n+e 1 ,l{x) 
\~ a n,rA-n+e r ,l{x) J 



Z = l,...,r. (2.9) 



From the lowering and raising equations stated in Theorems 12.11 and I2.2[ one readily 
derives differential equations of order r + 1 for the type II and the type I multiple orthogo- 
nal polynomials. As an illustration, we show how to derive differential equations satisfied 
by Pn{x) with r = 2, the idea of which can, of course, be extended to general situations. 



With r = 2 in (ESB, i.e., PJx 



(Pit(x), Pft-g^x), Ph-s 2 (x)) t , it follows from (12. 5ft and 



straightforward calculations that 



PUx) 

x) 



P" 



n 
pw 



X 



N{x)P n {x), 

N 1 {x)P n (x), 

N 2 (x)P H (x), 



where 



iVi(z) = N'(x) + N 2 (x), 

N 2 (x) = N"(x) + 2N'(x)N(x) + N(x)N'(x) + N 3 (x), 

and the derivative of a matrix valued function is understood in the entry-size manner. 
Hence, 



P' H {x) = (N) 00 (x)P n (x) + (N) 01 (x)P^ gl (x) + (iV) ( 
P%(x) = {N l ) 00 {x)P H {x) + (Ni) 01 {x)P n - gl (x) + (N 



P^'(x) = (N 2 ) 00 (x)P r 



)(a(x) Phs 2 (x), 
-Si [ x ) j 

\{x) + (N 2 ) 01 (x)P H - Sl (x) + (N 2 



Pfi— e 2 {X), 



(2.10) 
(2.11) 
(2.12) 



where we use the notation (M) ^ to denote the (i, j)-th entry of any given matrix M. Then, 
we can represent Pn-Sx an d Pn-s 2 m terms of Pft, P' H and P^ by solving (12.101) and (12. lip . 
Finally, replacing Pn-Sx and Pn-e 2 in (12.121) by these relations will lead us to the third 
order differential equations satisfied by P^. Since the exact formulas are cumbersome, we 
plan not to write them down here, but will present some concrete examples in Section [5j 



2.2 Compatibility conditions 

We finally state the compatibility conditions for the ladder equations for multiple orthog- 
onal polynomials. First we write the nearest-neighbor recurrence relations in a vector 
form. Let Ph{x) be defined by ( 12. 61) . then we have 

Pn+eAx) = W{n + e l] x)P H {x), (2.13) 
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for I = 1, . . . , r, where 



W(n + eu x) 



1 Bu(n) 







B,, 2 (n) 










'0"fi,r 






V i o 

is a matrix polynomial of degree 1 with 







(2.14) 



B hr (n)J 



The formula (12.131) follows from the nearest-neighbor recurrence relations (II. 6p . since by 
eliminating xP^(x) between Z-th and j-th relation in (11.61) . we have 

Pft+g^x) - Pn+e^x) = (bti,j - bfi t l)Pfi{x). 



Rewriting these relations and the Z-th relation in (ll.6p gives us (12.131) and (12. 14H . We then 
have 

Theorem 2.3. The compatibility conditions for the ladder equations stated in Theorem 
\2.1\ are given by 

N(n + ef,x)W(n + ef,x) = W'{n + ei,x) + Win + e x \ x)N(n; x), Z = 1, . . . ,r, (2.15) 

where N(n;x) is defined in (I2.5p . 

Clearly, by considering (I2.15P in an entry-size manner, one has at least r(r + l) 2 equal- 
ities, which are much more complicated than the single weight case. Similar compatibility 
conditions for ladder equations associated with the type I multiple orthogonal polynomials 
can be obtained easily by using (I2.15p . We omit the results here. 



2.3 Outline of the paper 

The rest of this paper is mainly devoted to the proofs of our theorems. We will prove 
Theorems 12.11 and 12.21 in Section [31 Unlike the treatment of the single weight case [TT] . 
Riemann-Hilbert (RH) problems and the Christoffel-Darboux formula for multiple orthog- 
onal polynomials will be two fundamental components in the derivations of the ladder 
equations; see also [27] for similar treatment to matrix orthogonal polynomials. We then 
derive the compatibility conditions in Section HJ where we also give a review of partial 
difference equations satisfied by the nearest-neighbor recurrence coefficients, obtained in 
[37] . for later use. We conclude this paper by applying our results to several concrete 
examples, namely, multiple Hermite and Laguerre polynomials, and multiple orthogonal 
polynomials with exponential weights and cubic potentials. 



3 Proofs of Theorems 2.1 and 2.2 



As mentioned before, the proofs of Theorems 12.11 and 12.21 rely on the Riemann-Hilbert 
(RH) problem and the Christoffel-Darboux formula for multiple orthogonal polynomials. 
In what follows, we first give a brief introduction of these two aspects. 



S 



3.1 Riemann-Hilbert problems and the Christoffel-Darboux for- 
mula for multiple orthogonal polynomials 

The usual orthogonal polynomials can be characterized by a RH problem of size 2x2 
[23]. It was shown by Van Assche et al. [39] that multiple orthogonal polynomials can 
also be described in terms of a RH problem, but for matrices of order (r + 1) x (r + 1). 

For convenience, we assume that all the measures are supported on an oriented, un- 
bounded and simple curve F in the complex plane. If the weight functions Wj are Holder 
continuous, we look for an (r+1) x (r + 1) matrix valued function Y satisfying the following 
RH problem: 

1. Y is analytic on C \ Y. 

2. Y possesses continuous boundary values Y + (x) (from the positive side of T) and 
YL(x) (from the negative side of T) , which satisfies 



Y+(x) = y_(ar) 



fl wi(x) w 2 (x) 

1 

1 

\0 



w r (x)\ 







1 



x g r. 
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3. As z — > oo, z G C \ T, we have 

Y(z) = {I + 0(1/ z)) diag^ 1 ™ 1 , z~ n \ z - n2 z' nr ). 



If the indices n, n — e*,-, j = 1, . . . , r, are normal, then there exists a unique solution in 
terms of the type II multiple orthogonal polynomials which is given by 



Y 



Pdz) 



where 



-2Tri-f Htr P H _ gr (z) 



1 



1 



dt 



-27ri<y HA P H _ Sl (z) -7^1 



2iri J T t — z 

Pn-g^wtit) 



t-z 



Pn-eMlU^t) 



t-z 



dt 



dt 



x n i- x p n - s m Wj {t) dt, 



~7rU 



Pn(t)w r (t) ^ \ 
P H „ Sl (t)w r (t) 



2ni ./r t — z 



t-z 



dt 



, Pn-g r (t)w r (t) 

-ln,r / 7 dt 

! T t-z 



j = l,...,r. 



/ 

(3.1) 

(3.2) 



There exists a similar RH problem for the type I multiple orthogonal polynomials to 
determine a matrix valued function X of dimension (r + 1) x (r + 1) such that 



1. X is analytic on C \ T. 
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2. For x G T, we have 



Xjx) =Xjx) 



( 1 

—Wi(x) 1 

— w 2 {x) 1 

^— w r (x) ••• 



... 0\ 

... o 


'•• 
1/ 



3. As z — > oo, z G C \ T, we have 



X(z) = (I + 0(l/z))di&g(z-W,z n \z n2 ...,z nr ). 

Assume that the indices n, n — ej, j = 1, . . . , r, are normal, then the unique solution X 
is given by 



X 



\2iri J T z — t 
where is defined in (jl.5p and 



gg(g 
^ - t 



dt 



2m 



z-t 



('n.r f Qn+e r (^) 



dt CfiiAfij r g li i[z) 



dt Cfi jT Afi + g r ^( K z) 



2mAff r (z) 



\ 



Cn,l-Afi+ei ,r \%) 



Cn,rAfi-{-er,r ( 



(3.3) 



(3.4) 



with Kfi+g-j being the leading coefficient of Aa+g.j. 

If the curve T is bounded, extra conditions at the endpoints of T are required in the 
RH problem to ensure the uniqueness of the solution, if the solution exists. 

A simple relation between the matrix functions X for the type I and Y for the type II 
multiple orthogonal polynomials is given by Mahler's relation (cf. [2E1 Theorem 23.8.3]): 



X(z) = Y~ T (z) 



(3.5) 



where A~ T is the transpose of the inverse of the matrix A. As we shall see later on, this 
relation plays an important role in our proof. 

We conclude this section with the Christoffel-Darboux formula for multiple orthogo- 
nal polynomials. This formula was given by Daems and Kuijlaars [17] as stated in the 
following theorem (see also [37] for another proof): 

Theorem 3.1. Suppose (fi$)i=o,i,...,[fi| a path in W such that no = 0, n\ft\ = n and for 
every i G {0, 1, . . . , |n| — 1} one has n i+ i — Hi = for some k G {1,2, . . . , r}. Then 



x 



\n\— 1 r 

V) P fti( X )Qn i+ l(y) = Pn(x)Qfi(y) - ^ a n,j P n-e,(x)Qn+e J (y)- (3.6) 

i=0 j=l 
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Observe that the right hand side is independent of the path (ni)i=o,i,...,\n\ i n A 
similar formula for multiple orthogonal polynomials of mixed type is derived in [TJ. This 
formula is of importance in the analysis of certain random matrices (TJ [30] and non- 
intersecting Brownian motions [18] . 

It is worthwhile to point out that if the system of measures is such that all the indices 
are normal (perfect system) and the functions 

{wi, xwi, . . . , x ni_1 Wi, w 2 , xw 2 , • • • , x n2 ~ 1 W2, • • • , w r , xw r , . . . , x nr ~ l w r } 

are linearly independent (this is exactly our assumption on the r measures in Theorem 
12. ip . then the Christoffel-Darboux formula also holds in a component- wise manner, i.e., 

\ri\-l r 

( X ~y)^2 P ni( X ) A n l+ i,k(y) = Pn(x)A Hjk (y) - ^ ^jPn-e, (x)A H+gj >k (y) , (3.7) 
i=0 j=l 

for y E supp(/ifc) and k = 1, . . . , r. To see this, we rewrite ( 13. 6ft as 

r 

*Y^R k (x } y)w k (y) = 0, 

k=l 

where 

\H\-1 r 

R k {x,y) = (x-y)^2 p nX x ) A n t+1 ,k{y) - Pn{x)A^ k {y) + ^2an,jPn-s j {x)A n+Sj , k {y). 

i=0 j=l 

With x fixed, it is easily seen that each Ri(x,y) is a polynomial in y of degree less or 
equal to n k . The linear independence then implies Ri(x,y) = 0, which is (13. 7p . 

3.2 Proof of Theorem 1231 

Let us define an (r + 1) x (r + 1) matrix M by 

M(x) = M{n; x) = (M ij (x)) < ij < r := Y'{x)Y-\x), (3.8) 

or, equivalent ly, 

Y\x) = M(x)Y(x), (3.9) 

where Y is given by (13. ip . Since the first column of Y is given in terms of the type II 
multiple orthogonal polynomials, it is easily seen that 

r 

Pt(x) = M 00 (x)P H (x) - 2m-y itJ Mo j (x)P it .g. (x), (3.10) 

j'=i 



and 



Pk-sM = -W^— \M i0 (x)P n (x) -^m-ynjMijWPrt-zix)) , (3.11) 



for i — 1, . . . , r 
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In the following we use the notation 

C(f)(x)= [p£-dt, xeC\T, 
Jy t — X 

for the Cauchy transform of /. The derivative dC(f)(x)/dx is denoted by C(f)'(x). 
On the other hand, the Mahler relation (13. 5p gives 

M(x) = Y'(x)X T (x). 

It then follows from (13. ip and (13.31) that 

r 

M 00 (x) = -P' R {x)C{Q R ){x) + ArtAx)C(PnW k )'(x), (3.12) 

k=l 

M 0j {x) = g (-PL(x)C(Q n+g .)(x) +^A a+gjik {x)C(P H w k y(x)\ (3.13) 
with j = 1, . . . , r. For i, j — 1, . . . , r, we have 

M i0 {x) = 2m lHti ^P' n _ Si {x)C{Q n ){x) - ^A H , k {x)C{P n _ Si w k )'{x) j , (3.14) 

My (a;) = ln,iCn,j ^P n _g i (x)C(Qft + g j )(x) - y~) A i t + g. }k (x)C(Pn-g i w k )\x)\ . (3.15) 

The entries in M can be simplified in view of the following two elementary facts. 
Proposition 3.1. If 

J t h f(t)dt = 0, k = 0,...,n — l, 

then we have 

p(x)C(f)(x) = C(pf)(x) 
for any polynomial p of degree less than or equal to n. 



Indeed, 



and (p(t) — p(x))/(t — x) is a polynomial of degree less than or equal to n — 1, hence, by 
assumption, the integral disappears. 

Integration by parts gives us the following statement. 

Proposition 3.2. If f is a differentiable and integrable function that vanishes at the 
endpoints ofT, one has 

C(f)'(x) = C(f')(x). 
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Hence, the orthogonality of P and Q (see (jl.3p and (jl.4j) ) implies 

Pi(x)C(Q n )(x) = C(PLQ H )(x), (3.16) 
A^ k {x)C{P H w k )'{x) = C(Art, k P n w k )'(x) - C(P H w k )(x)A' H!k (x), k = 1, . . . ,r. (3.17) 

Inserting the above formulas into (13.121) . it is readily seen that 



Moo(x) = -C{P' H Q H ){x) + {C(A Htk P n w k )'(x) - C(P H w k )(x)A' Hjk (x)) 

r 

= -C{PLQn){x) + C(P H Q H )'(x) - £ C(P^ fc )(x)A'- fc (x) 

r 

= C(P^)Or) - ^C(P^ fe Ay (: 
fc=i 

= c(p h J2 A n,kw' k ) (x) = -C [ P n A n,kV k w k ) (x) , (3.18) 



k=l 
(x) 



k=l / \ k=l 

where v k (x) := —lnw k (x). In a similar way, we have 



M i0 (x) = 27Tij n ,iC [Prt- gi A n,kV' k W k ) (x). (3.20) 



k=l 



The simplification of My with i, j > 1, however, is a little bit tricky. The problem 
is that Aft +gjik (x)C(Pfi-g i w k )'(x) in (I3.15P can not be written in a form similar to (I3.17P 
(indeed, for i = j = k by comparing the degrees of A^ + g jtk and orthogonality relations 
( II. 3p . we see that assumptions of the first proposition are not fulfilled). To overcome this 
difficulty, recall that YY^ 1 = YX T = I, the entry on both sides gives 



7: 



n,iCri,j Pn-eA X ) C (Qfi+? 3 )( X ) ~ ^ A H+ g^ k {x)C (Pfj-e^k) (x) = 5ij, l,j > 1. 



k=l 



Note that P^g i C(Qfi + g j ) = C(P f i^g i Qfi + g j ). Taking the derivative of the above formula 
with respect to x, we see 



A H+ g^ k (x)C{P^g l w k )'{x) = Y A ft+ g, )k (x)C(Pft-s i w k )(x) 

k=l k=l 

~ & (Pri-eiQn+eMx) — C '(Pfi^g.Q ft + g.) (x) 



(3.21) 



Using ( KTDf in flBTTo} we have 

r 
k=l 

= ^fn,iCri,jC ^ Pjt-gjAfi+g.fiV'kWk (x), i,j > 1. (3.22) 



v jfc=l 
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Replacing M tJ in (I3TTUD and (I3~TTD by f[3~TgD - (l3T2TIj) and (13321) . it follows 
P' H {x) = -C [PnJ2A a>k v' k w k ) (x)P n (x) 



k=l 



(3.23) 



and 



(3.24) 



+ ^2 a n,jC I Ptiy^ J Afj + g hk v k w k J (x)P f i-s j (x), 
j=i V k=l J 

P'n-e^) = -C (^n-e^A^ k {x)v' k W k ^l (x)P a (x) 
r / r \ 

+ QfijC I Pfj-Sj An+s hk v' k w k J (a^P^-e^x), 

j=l \k=l ) 

for i = 1, . . . , r, where we have used the fact that 

Cnj^/nj = ®"ri,ji (3.25) 

see (PUT), (GOD, (HUDD and f l3^|) . 

To show (12.11) . we observe from (13.231) that 



PL(x) = P n (x) [ Pn(t)J2A H 4t) 
Jr k=i 



v' k {t)-v' k {x) 



x — t 



w k (t) dt 

X(*)-^(z). 



- anjPfi-g^x) [ Pn(t) ^ A H+gj , k (t) " kK "' w k (t) dt + E(x), (3.26) 

3=1 Jr k=i x 1 



where 



E i x ) = ^2 [ ^2 a n,jC{P H A ri+S: . tk w k ){x)P H ^ :j {x) - C(P H A Htk w k )(x)P r -i(x) J v' k (x). 
k=i \j=i J 

(3.27) 

It is then equivalent to show E = 0. To this end, we apply the Christoffel-Darboux 
formula (13.71) to each coefficient of v' k in (I3.27P and obtain 



^2a H jC(P f tA H+ g j:k w k )(x)P H _g j (x) - C(P H Afi >k w k )(x)P a (x) 
3=1 



(y ' Pn{t)A Hi+uk {t)w k {t)dtjP n 



where we have used the orthogonality of Pfi (I1.3P in the last equality. Hence, it is imme- 
diate that E = and Q2HD follows. 

In a similar way, we see from (13.241) that 

P' a _ Si {x) = P H (x) [ P H ^{t)j2An,k{t) V ' kit) - V ' k(x) w k (t) dt 
JL k=l 

T««/^W / Pn-dt) Y^A H+ ^ k {t) v ' k{t) - v ' k{x) w k (t) dt 

3=1 J ? k=l X t 

+ E i {x), (3.28) 
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where 



E i( x ) = ^2 [ ^2 a rtdC{P ri _g i Af i+e ,^ k w k ){x)P H ^g j {x) - C(P a ^g i Afi tk w k )(x)P n (x) ) v' k (x). 
k=l \j=l J 

Again with the aid of the Christoffel-Darboux formula (13.71) and the orthogonality of 
Pn-Sii it is readily checked that 



^ ] a fj,jC'{Pfl-S i ^n+ej,k w k) { x )Pfj-gj j x ) ~ C (Pfi-g i A^j,W k )(x)Pfi(x) 



^ ^ Ps-e-(*)Aj i+ i,fe(^)w fe (t)cf^ -P«i(a:) = Ph-sA^^j- 



This, together with (13.261) . implies 



(3.29) 



Combining (13.281) and (I3.29j) then leads to (12. 2p . This completes the proof of Theorem 

o 

Remark. There is a simple relation between the matrices M in (13.81) and iV in ( 12. 3ft . 
Note that the first column on both sides of (13.91) gives us 



/ Pn( 



d_ 

dx 



x 



\ 



-2m^iP H _ Sl (x) 



M(x) 



\-2-ni~f H)T P H _ Sr ( x ) j 
We then obtain from (12. 5ft that 

1 1 



/ Pn{x) \ 
-2mj H;1 P n ^ gl (x) 

\-2-Ki^fi,rPn-S T { x ) J 



N(x) = diag 1 



2m-in,i 



M(x) diag (1, -2^7^,1, • • ■ , -2-K%^ n . 



(3.30) 



3.3 Proof of Theorem ¥Z72\ 

We shall prove Theorem 12.21 by establishing (I2.9p . The formulas (12 ,7p and (12.81) then 
follow from straightforward calculations using (12. 9p and (12. 3p . 

To see (I2.9p . we first take a derivative with respect to x on both sides of the equation 
Y{x)Y~ l (x) = I and obtain from (13 .8p that 

Y(x)(Y- 1 Y(x) = -M(x). 

Equivalently, in view of (13. 5p . 

(X T ) f (x) = -X T (x)M(x). (3.31) 
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Next, we observe from (I3.3P that, by transposing the l-th (I = 1, • • • , r) row on both sides 
of fl3T3TD . we get 



d 

dx 



( 2mA^i{x) \ 

Cn,1^4n+e*i,i(- r ) 
\ c fi,rAn+e r ,l{ x ) j 



-M T (x) 



( 2mA^(x) \ 

Cn,1^4n+e*i,«( a ') 

\Cfi )r Art + g r j(x) J 



that is, 
d 

dx 



( A nA X ) \ 

\A R+ g r j(x)J 



diag ( — , — , • • • , — ) M T (x) diag (27rz, ch,i, • • • , c^ r ) 

4 ^TTi Cn,l Cn,r , 



X 



\An+e r ,l( X )J 

Replacing M T in the above formula with the aid of (I3.30p . we finally arrive at (12 .9p after 
simple manipulations using (13. 25ft . This completes the proof of Theorem 12.21 



4 Proof of Theorem 2.3 



To show (I2.15p . we note that equations (I2.5P and ( 12 . 131) can be viewed as a Lax pair. 
Differentiating both sides of ( 12 . 1 3j) . together with (I2.5p . implies 

N(n + ei;x)P H+gl (x) = (W'(n + ef,x) + W(n + e l ;x)N(n;x))P H (x), 

for I = 1, . . . , r. Using (I2.13P again, it follows 

N(n + eu x)W{n + e,; x)P H (x) = (W'(n + e { ; x) + W(n + e x \ x)N(n; x)) P s (x). (4.1) 

It is now sufficient to show that the components of the vector are linearly independent 
whenever n is a normal index. This assertion is already shown in [37] . For the convenience 
of the reader, we repeat the argument here. Suppose that (cq, ci, . . . , c r ) are such that 

r 

CqPti + CjPn-ej = 0, 

i=i 

then by comparing the leading coefficients, it follows cq = 0. If we multiply by x nk ~ 1 and 
integrate with respect to fi k , then 



c k J x nk x Pfi-g k {x) dfi k (x) = 0. 



We claim that the above integral does not vanish. Indeed, if the integral vanishes, then 
Pfi — aPa_g k is a monic polynomial of degree \n\ satisfying the orthogonality relations 
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(jl.3p for any a, which contradicts the normality of n. As a consequence, = for 
A; = l,...,r. The linear independence of the polynomials in and (14.11) gives the 
compatibility conditions (I2.15p . This completes the proof of Theorem 12.31 

Remark. From (I2.13p . it is readily seen that one can obtain the vector P^ + g. + g. in two 
different ways. The first way is to first compute Pn+e; from P^ and then to compute 

Pn+ei+e, from Pfi + g i : 

P ' n+Si+ g.(x) = Win + e*i + e}; x)W(n + x)Pn(x). 

The second way is to first compute Pa+gj from P^ and then to compute P^ + g i+ g, from 
- ■ 

1 n+ej ■ 

P H+ g. + g j (x) = W(n + ej + ti\ x)W(n + ef x)P n (x). 
Hence, we obtain another kind of compatibility condition 

W{n + e*j + tj\ x)W(n + ef, x) = Win + e} + ef, x)W(n + e}; x). 

This was observed in [37], and from the explicit formula of W the following nearest- 
neighbor recurrence relations for multiple orthogonal polynomials were obtained. 

Theorem 4.1. [37, Theorem 3.2] Suppose all multi-indices n e N r are normal. Suppose 
1 < i 7^ j < r , then the recurrence coefficients for the nearest neighbor recurrence relations 
(CQ|) satisfy 



k=l 





bfL-\-£j,i bfi,ii 


r 

~ ^ ] a ft+ei,k 
k=l 


\ u n+Si,j u n,j/ 




bft—ei,j bfi—ei,i 


®ri+ej ,i 





This theorem tells us that the recurrence coefficients of multiple orthogonal polyno- 
mials cannot be arbitrary, but satisfy a system of partial difference equations. This is not 
the case for the usual orthogonal polynomials. 



5 Some examples 

In this section, we shall derive the ladder equations and differential equations for several 
examples of multiple orthogonal polynomials. For convenience, it is assumed that r = 2 
throughout this section. The (normal) multi-index n is now given by (n, m) G N 2 . We 
also use the following notation for the recurrence coefficients: 

The recurrence relations (11. 6p then read 

3jPn,m\X) Pn,m+1 \p^) d n ^ m P n ^ m iyX} -\- O nm P n —\ m ix^ -\- b n ^ m P nrri —\ix^ (5.2) 
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with ao, m = and b nt o = for all n, m > 0. In view of Theorem 14. 1[ the following relations 
hold: 

(5.3) 
(5.4) 

(5.5) 

(5.6) 

^n,m—l u n,ra— 1 

5.1 Multiple Hermite polynomials 

Multiple Hermite polynomials H n m are type II multiple orthogonal polynomials defined 
by 





7 

Q"n,m 




"■n+l,ra 




(j^g^ j dn+l,m ^n,m j 

\Cn,m+l Cj^m / 






Cn,m d n ^ m 




bn+l,m 


r A ' 
<^n— l,m "71—1,7)1 

Cn,m ^n,m 




bn,m 


Cn,m— 1 d nm —\ 



oo 

71,77) 

-OO 
OO 



x k H nm (x)e x2+cix dx = 0, A; = 0,1 ra - 1, 



x k H n>m {x)e~ x +C2X cfe = 0, fc = 0, 1, . . . , m - 1, 

where c\ ^ c 2 ; cf. [8], [25J §23.5] and [3S1 §3.4]. Hence, the functions Vk(x) in Theorem 
12. II are equal to x 2 — CkX. A simple calculation with the aid of ( 11.3jl then gives 




iV(n, m; ac) = -2 2x - c x . (5.7) 



Recall that N(n,m;x) is defined by (I2.5P and explicitly given in (j2.4p . This, together 
with ( 12.141) and the compatibility conditions ( 12.151) . leads to the the following system of 
difference equations for the recurrence coefficients: 

(5.8) 
(5.9) 
(5.10) 
(5.11) 
(5.12) 
(5.13) 
(5.14) 

It turns out one can easily obtain the recurrence coefficients explicitly from these equa- 
tions. Indeed, from (15.101) . it follows 

Cn,m Ci/2, d nm C2/2. 

Since c\ 7^ c 2 , equations (15.111) and (15.121) imply that b n ^ m is independent of n, while (15.131) 
and (15. 14j) show that a n>m is independent of m. Hence, by (15.81) and (15. 9p . it follows 

2a n ,m — 2a n)m +i + 1 = 0, 2b njTn — 2b n +\ )m + 1=0. (5.15) 



2dn.ra 2d n 


+ 1,77) + 


2b nrn 2b n +\ m + 1 


= 0, 


2dn,m 2a n 


,7)1+1 + 


2b nm 2b nm +i + 1 


= 0, 






2c njm Ci, 2d nm 


= c 2 


2(d nm . 


-1 _ C 7) 


,771—1) (02 2c nim ) 


= 0, 


2&n+l,77) {dn, m— 1 


£-71,771—1 


) b nm (c2 20^^) 


= 0, 


2(4- 


-1,77) ^Tl— 1,77)) + Cl 2d nm 


= 0, 


2*3-7), m+l (dn— 1,771 


Cn— 1,77), 


) + ^ri,m (C-l 2d n m ^) 


= 0. 
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Note that our initial conditions are a 0jm = and b Ut Q = 0, we then obtain from (15.151) 

Thus, we arrive at the following lowering and raising equations for multiple Hermite 
polynomials /f n , m : 

), (5.16) 

H n-i,m( x ) = -2H n>m (x) + (2x - ci)# n _i im (x), (5.17) 

H' n , m -M = ~2H n , m (x) + (2x - c 2 )H n>m ^{x). (5.18) 

Oberve that equation (15.171) (similarly (15.181) ) can be alternatively written as 



(e- x2+c ^H n _ hm (x))' = -2e- x2+c ^H n ^{x). 



These formulas are not new but can already be found in [28, §23.8.2]. 

Finally, by (1 5. 7ft and arguments at the end of Section 12.11 we see that the type II 
multiple orthogonal polynomial H n m satisfies the following linear differential equation of 
order 3: 

p"'(x) + (ci + c 2 - Ax)p"{x) + (ci(c 2 - 2x) + 2(m + n - 1 - c 2 x + 2x 2 )) p'(x) 

+ 2{c\m + c 2 n — 2(m + n)x)p(x) = 0. 

Similarly, we can derive a third order differential equation for the type I multiple Hermite 
polynomials A( n>m \i, I = 1,2, which is given by 

q"'(x) - (ci + c 2 - Ax)q"(x) + (ci(c 2 - 2x) + 2(m + n - 1 - c 2 x + 2x 2 )) q'(x) 

— 2{c\m + c 2 n — 2(m + n)x)q{x) = 0. 

We point out that the above differential equation is independent of /, which can also be 
seen from (12. 9j) . 

5.2 Multiple Laguerre polynomials of the second kind 

These polynomials are defined by the orthogonality conditions 

x k L n , m (x)x a e~ ClX dx = 0, k = 0, 1, . . . , n - 1, 
x k L nm (x)x a e~ C2X dx = 0, k = 0, 1, . . . , m - 1, 



o 







where we assume that a > and ci,c 2 > with c\ ^ c 2 ; cf. [8], [331 Remark 5 on 
p. 160], [2H1 §23.4.2] and [38], §3.3]. Hence, the functions Vk(x) in Theorem 12.11 are equal 
to —a Inx + CfcX, k = 1,2. An appeal to the biorthogonality conditions (11.41) implies that 



N(n, m;x) = — 
x 



( -a f™ Ln > m{t)< ? n ' m{t) dt A(n,m) B(n,m) \ 

-a J °° L -vnfflQn, m (f) dt C (n,m) -a + Cl x D{n,m) 

v -ajr o °° L "'"- l(t t )Q "' m(t) dt E(n,m) F(n,m) - a + c 2 x j 



(5.19) 
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where Q n ,m is defined in (jl.5p with Wj(x) 



x a e CjX , j = 1,2, and 



A(n, m 
B(n, m 
C(n, m 
D(n, m 
E(n, m 
F(n, m 



aa r 



ab r , 



aa r 



ab r . 



aa r 



ab r , 



(t)Qn+l,m{t) 
t 

-^n,m (t)Qn,m-\-l \t) 



dt. 



dt. 



J n—l,m 



(t)Qn+l,m\t) 



J n—l,m 



(t)Qn,m+l(t) 



Ln,m-l{t)Qn+l,m{t) 



J n,m—1 



{t)Qn,m+l{t) 



dt, 



dt, 



dt, 



dt, 



(5.20) 
(5.21) 
(5.22) 
(5.23) 
(5.24) 
(5.25) 



are certain constants depending on n and m. We want to give an explicit representation 
of N(n,m;x) in (15.191) . We first observe that the first column in (I5.19P can be evaluated 
by comparing the leading coefficients on both sides of the differential equations (12. ip and 
(P . For instance, by and (15TT9|) . it follows 



-a 



J -'n,rn {t)Qn,m{t) 



dt I L r 



(x 



+ A(n, m)L n _ ljm (x) + B(n, m)L n>m _i(x). 

Recall that L nm (x) = x n+m + . . ., hence by comparing the coefficients of order n + rnon 
both sides of the above equation, it is readily seen that 

J -'n,m (t)Q n ,m(t) 



a 



-dt 



Similarly, we obtain 



a 



dt — ci, a 



oo L 



(n + m). 

n,m— l{t)Qn,m{t) 



t 



dt = c 2 . 



(5.26) 



(5.27) 



To estimate the integrals in (15. 201) — ( T5.25j) . we shall make use of the compatibility condi 
tions (12.151) and the known results of the recurrence coefficients: 

(n + m + a)n (n + m + a)m 



and 



2n + m + a + l m 

1 , d r< 

ci c 2 



r ,2 
c 2 



n + 2m + a + 1 n 

+ -; 

C 2 Ci 



(5.28) 
(5.29) 



see [371 §5-4]. The compatibility conditions ( 12 . 1 5[) in this case are given by 

-K 

xN{n + 1, m; x) 




I X 

















H 


s 











xN(n,m;x), (5.30) 



1 Cn ,m—l, 
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and 



d n ^ m Q"n,m bn,m 



f X 





°\ 













H 


\o 










xN(n,m + 1; a;) | 1 c n _i, m - d n _i iTn 

1 

1 Cn—l,m ^n— l,ra 

I xN(n,m;x). (5.31) 
1 

With the aid of (I5.19P and the explicit formulas (l5.26jHI5.29p . we see from the (2, 2)-entry 
in flB^Ojl that 

(n + m + a)n 
A(n,m) = cxa^m = . 

Cl 

Similarly, the (3, 3)-entry in (I5.3ip implies 

% , (n + m + a)m 

B(n,m) = c 2 b n , m = . 

c 2 

Then, we observe from the (2, 3)-entry in (I5.30p that 

cA,« + (rf n ,m-i - c n , m _i)-D(n + 1, m) = B(n, m), 

which gives 

JJ{n } m) = . 

c 2 

Similarly, the (3, 2)-entry in (I5.3ip implies that 

h{n,m) = , 

Cl 

which is independent of m. Finally, we obtain from the (2, l)-entry in (I5.30P that 
c\Cn,m + C(n + 1, m) — a + D{n + 1, m) = n + m, 

hence, 

C(n, m) = —n. 
An appeal to the (3, l)-entry in (15.301) gives 

c 2 c n , TO + E(n + 1, m) + F(n + l,m) - a = n + m + c 2 (c„ im _i - d„ jm _i), 

thus, 

F(n, m) = —m. 

Combining all these results, we find the following lowering and raising equations for mul- 
tiple Laguerre polynomials of the second kind L n , TO : 

ti , \ / \ T i \ f(n + m + a)n\ , 
xL nm [x) = (n + m)L n>m (x) + I 1 L n _ lym {x) 

( (n + m + a)m\ 
+ ( " —j L^m-^x), 

xL 1 m (x) = —c\L nm (x) — (n + a — c\x) L„__i m (x) L nm _i(x), 

c 2 

a; -^nm-i( ;r ) = — c 2 L nm (x) L n _i m (x) — (m + a — c 2 x) L nm _i(x). 

Ci 
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Furthermore, the third order differential equation satisfied by L n>m is given by 

x 2 p"'{x) - (x 2 (c! + c 2 ) - 2x{a + 1)) p"(x) 

+ [x 2 c x c 2 - x[(ci + c 2 )(a + 1) - nci - mc 2 ] + a(a + 1)) p'(x) 
— [xc\C 2 {n + m) — a{nc\ + mc 2 )) p(x) = 0, 

which can also be found in [51 §4.3]. Similarly, it can be shown that the type I multiple 
Laguerre polynomials of the second kind satisfy the following differential equation: 

x 2 q"'(x) + (x 2 (d + c 2 ) - 2x(a - 1)) q"(x) 

+ {x 2 C\C 2 — x[(ci + c 2 )(a — 1) — nci — mc 2 ] + a(a — 1)) q'(x) 
+ (xcic 2 (n + m) — a{nci + mc 2 )) q(x) = 0. 



5.3 Multiple Laguerre polynomials of the first kind 

These polynomials are defined by the orthogonality conditions 



x L nm (x)x ai e x dx = 0, k = 0, 1, . . . , n — 1, 



x L n ^ m (x)x a2 e x dx = 0, k = 0, 1, . . . , m — 1, 



where ai,a 2 > and a\ — a 2 ^ Z; cf. [8], [281 §23.4.1] and [381 §3-2]. The functions Vk{x) 
in Theorem 12.11 are now equal to — c^lnx + x, k = 1,2. A straightforward calculation 
using ( 11. 4fl gives 



iV(n, m;x) = — 
x 



where we define 



^ _ joc L„, m (t)Q», m (t) ^ A , m) 

f °° ^-i.mWQn,m(«) dt C*(n,m)-a 1 + x 



Wo' 



OO L n ,m-l(*)Qn,mffl 



E*(n, m) 



B*(n,m) \ 
D*(n,m) 
F*(n, m) — a 2 + x J 



Qn,m(*) = ("lAn,-),lW tQl + a 2 A (n ^ 2 (t)t a2 ) e~\ 

with y4( nm )j, j = 1,2 the associated type I multiple orthogonal polynomials and 

Ln,m\t)Q n -i-i <m {t) 



A*(n, m 
B*(n, m 
C*(n, m 
D*(n, m 
E*(n, m 
F*(n, m 



t 



dt, 



L n ,m [t) Qn,m+1 W 



dt, 



Ln-l,m{t)Qn+l,m(' t ) 
t 

Ln-l,m(t)Qn tm+ i(t) 
t 

Ln,m-l{t)Q* n+lm {t) 
t 

Ln,m-l(t)Qn jm +l(t) 



dt. 



dt. 



dt. 



dt, 
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are certain constants depending on n and m. Again by comparing the leading coefficients 
on both sides of the differential equations (12. ip and (I2.2p . we have 

J -'n,m 

at = —{n + m) 

^n-l,m(Wn, mW ,^ _ f°° ^ro-lW^mW dt ~ 1 



t 

Since the recurrence coefficients are explicitly given by 

_ n(n + oti)(n + ati — a 2 ) _ m(m + a 2 )(m + a 2 — «i) 

n — m + oti — a 2 m — n + a 2 — aci 

c n ,m = 2n + m + «! + 1, d njm = n + 2m + a 2 + f , 

see [371 §5.3], the same strategy as in Section I5T21 gives 

m) = n(n + 0^ + 0,-0^ ^ m) m(m + a 2 )(m + a 2 - a x ) 



n — m + ai — a 2 m — n + a 2 — oc\ 

C*(n,m) = —n, F*(n,m) = —in, 

D*(n,m) = E*(n,m) = 0. 

Hence, we find the following ladder equations for multiple Laguerre polynomials of the 
first kind L njTn : 

, / n(n + aci)(n + a% — a 2 ) . , 

xL nm {x) = [n + m)L n<m {x) + ( ■ ) L n _ hm (x) 

n — m + a.\ — a 2 



Ln,m—l \pC)i 



m(m + a 2 )(m + a 2 — ai) 
m — n + a 2 — ot\ 
xL n _ lm {x) = —L n ^ m {x) — (n + oil ~ x) L n _\^ m (x), 
xL' nm _i(x) = -L n>m (x) - (m + a 2 -x) L„ jm _i(x). 

Furthermore, the third order differential equation satisfied by L n ^ m is given by 

x 2 p"'(x) + (-2x 2 + (ai + a 2 + 3)x) p"(x) 

+ (x 2 — x(ai + a 2 — n — m + 3) + (ai + l)(a 2 + 1)) p'(x) 
— (x(n + m) — (n + m + nm + aim + a 2 n)) p(x) = 0, 

which can also be found in [5, §4.3]. Similarly, we have that the type I multiple Laguerre 
polynomials of the first kind satisfy the following differential equation: 

x 2 q"'{x) + (2x 2 - (ai + a 2 - 3)x) q"(x) 

+ (x 2 — x(cti + a 2 — n — m — 3) + (ai — 1)(«2 — 1)) q'{x) 
+ (x(n + m) — (mn — n — m + aim + a 2 n)) q(x) = 0. 
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5.4 Multiple exponential polynomials with cubic potentials 

These polynomials are defined by the orthogonality conditions 

J x k P n , m (x)e- x3/3 - ClX dx = 0, k = 0, 1, . . . , n - 1, 
x k P nm (x)e- x3/3 - C2X dx = 0, k = 0,1 m - 1, 



5 

r 



where c\ ^ c 2 and the contour F is taken to be the set {x : 9ft x 3 < 0}, or simply, 
{x : argx = ±27r/3}. This definition generalizes the orthogonal polynomials introduced 
by Magnus in [32} §6]. Since Vk(x) = x 3 /3 + CkX, we obtain from (11. 4ft and (ll.6p (or (15.1 ft 
and ( 15. 2 p in case r = 2) that 



2 



iV(n, m; x) = ( -x - c n _i, m x 2 + a„ jm + c a 6 n , m 

If we introduce the matrix form for N: 

N(n, m; x) = F 2 (n, m)x 2 + F\(n, m)x + Fq(ti, m) 
and similarly for W in ( 12.140 : 

W(n + 1, m; x) = i?x + W(n, m + 1; x) = i?x + -R2, 



where 



'1 X 

/? : 10 

,0 0, 



and 





( Cn,m 




Ri = j 















J , Ri = I 1 c n _i i?Tl — rf n _ 1]?n 

n,m— 1 Cn,m— 1/ V ^ 

then the compatibility conditions ( 12 . 1 5[) are given by the following equations (which are 
not identically zero), after collecting the coefficients with respect to x: 

F x {n + 1, m)Ri + F (n + 1, m)R = RF (n, m) + RxFi{n, m), 
Fi(n, m + 1)R 2 + F (n, m + l)R = RF (n } m) + R 2 Fi(n, m), 

F (n + 1, m)Ri = R + RiF (n, m), 

F {n, m + l)R 2 = R + R 2 F (n, m). 

Although in this case we do not know the exact expressions of the recurrence coefficients, 
the above relations indeed give us many nonzero equations for these coefficients. It can 
be shown that all but the following four equations can be simplified using them and ( 15. 3p 
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and fl£H) : 




(5.32) 
(5.33) 



(5.34) 



(5.35) 



In addition to equations fl5.3p - fl5.6p one can get, for instance, 



det 




In particular, we see that equations (15.321) and (I5.33P are nonlinear, which is different 
from the previous examples. We hope these relations will be helpful in the further study 
of multiple exponential polynomials with cubic potentials. 

Finally, we point out that it is also possible to derive the differential equations for the 
associated type I and type II multiple orthogonal polynomials. Since they are cubersome, 
we shall not write them down here, but mention that they are of the form 



where dk{x) are polynomials in x of degree k with the coefficients depending on the 
recurrence coefficients a n<m , b n , m , c n<m , d n>rn . 
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